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Abstract: We perform a detailed dynamical analysis of generalized Galileon cosmology, 
incorporating also the requirements of ghost and instabilities absence. We find that there 
are not any new stable late-time solutions apart from those of standard quintessence. 
Furthermore, depending on the model parameters the Galileons may survive at late times 
or they may completely disappear by the dynamics, however the corresponding observables 
are always independent of the Galileon terms, determined only by the usual action terms. 
Thus, although the Galileons can play an important role at inflationary or at recent times, 
in the future, when the universe will asymptotically reach its stable state, they will not 
have any effect on its evolution. 
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1 Introduction 

There are two approaches one can follow in order to describe the observed universe accel- 
eration. The first is to introduce the concept of dark energy, usually adding extra scalar 
fields (see [1] and references therein) in the right-hand-side of the field equations of Gen- 
eral Relativity, while the second is to modify the left-hand-side of the general relativistic 
field equations, that is to modify the gravitational theory itself (see [2, 3] and references 
therein) . 

A recently re-discovered, very general class of scalar-field theories is based on the in- 
troduction of higher derivatives in the action, with the requirement of maintaining the 
equations of motion second order. Although the most general second-order theories avoid- 
ing the Ostrogradsky instabilities [4] were already derived in [5], a particular class, dubbed 
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Galileon, was constructed in [6] for the Minkowski metric and in [7, 8] for dynamical geome- 
tries, while in [9] the results of [5] were re-discovered in the context of (extended) Galileon 
framework. In this formalism the Lagrangian is suitably constructed in order for the field 
equations to be invariant under the Galilean symmetry <f) — > <j> + c, d^(j) — > d^(j) + in the 
limit of Minkowski spacetime, with c, b^ constants. The four-dimensional Lagrangian that 
preserves these symmetries contains five unique terms consisting of scalar combinations of 
6^(1), d^d u 4> and (since General Relativity does not accept the satisfaction of Galilean 
symmetry in curved spacetime, higher derivatives are necessary in the action [10-12]), and 
the corresponding couplings are not renormalized by loop corrections [13]. Furthermore, a 
significant advantage is that the derivative self-couplings of the scalar screen the deviations 
from General Relativity at high gradient regions (small scales or high densities) through 
the Vainshtein mechanism [14], thus satisfying solar system and early universe constraints 
[15-17]. 

Application of the above construction in cosmological frameworks gives rise to the 
Galileon cosmology, which proves to be very interesting and has been investigated in detail 
in the literature. In particular, one can investigate the universe evolution and late-time 
acceleration [18-27], inflation [28-35] and non-Gaussianities [36-40], the reheating of the 
post-inflationary universe [41], the growth history [42-44], the cosmological bounce [45, 46], 
the cosmological perturbations [47-51], the spherical solutions [52-57], the stability issues 
[58-63] , and also he can use observational data in order to constrain the parameters of the 
theory [64-72]. Moreover, one can study and extend the properties of Galileon theory itself 
[73-93], and examine the relation of Galileons with other frameworks [94-105]. 

Since Galileon cosmology exhibits interesting phenomenological features, in the present 
work we perform a phase-space and stability analysis of such a scenario, investigating in 
a systematic way the possible cosmological behaviors, focusing on the late-time stable 
solutions. Such an approach allows us to bypass the high non-linearities of the cosmolog- 
ical equations, which prevent any complete analytical treatment, obtaining a qualitative 
description of the global dynamics of these models, which is independent of the initial 
conditions and the specific evolution of the universe. Moreover, in these asymptotic so- 
lutions we calculate various observable quantities, such as the dark energy density and 
equation-of-state parameters and the deceleration parameter, and in order to ensure that 
these solutions are free of ghosts and Laplacian instabilities we additionally calculate the 
relevant perturbations quantities. Interestingly enough, our analysis shows that Galileon 
cosmology does not exhibit any new stable late-time solutions apart from those of standard 
quintessence, and moreover the corresponding observables are always independent of the 
Galileon terms. Thus, although the Galileons can play an important role at inflationary or 
at recent times, in the future, when the universe will asymptotically reach its stable state, 
they will not have any effect on its evolution. 

The plan of the work is the following: In section 2 we briefly review the Galileon 
cosmological paradigm and in section 3 we perform a detailed phase-space analysis. In 
section 4 we discuss the cosmological implications and the physical behavior of the scenario. 
Finally, in section 5 we summarize the obtained results. 
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2 Generalized Galileon cosmology 



In this section we briefly review Galileon cosmology for the most generalized scenario, 
presenting the background cosmological equations and the conditions for the absence of 
instabilities [9, 18, 77]. As it is known, in order to avoid the Ostrogradski instability [4] 
it is desirable to keep the equations of motion at second order in derivatives, and thus the 
most general 4-dimensional scalar-tensor theories having second-order field equations are 
described by the Lagrangian [18] 

5 

C = Y J ^, (2-1) 

i=2 

where 

C 2 = K(4>,X), (2.2) 

£ 3 = -G 3 (0,X)n0, (2.3) 

U = G 4 (0, X) R + G A , X [(D0) 2 - (V„V^) (V"V»] , (2.4) 
C 5 = G 5 (<P,X)G I1U (V^V» 



- i G 5 ,x [(D0) 3 - 3(D0) (V„ V„0) (V"V» + 2(V^V a 0) (V Q V^) (V^V A 



(2.5) 



The functions K and Gj (i = 3, 4, 5) depend on the scalar field <j) and its kinetic energy 
X = —d^4>d l j i (j)/2, while R is the Ricci scalar, and G^ v is the Einstein tensor. Gi t x and G^ 
(z = 3, 4, 5) respectively correspond to the partial derivatives of Gj with respect to X and 
cp, namely G^x = dGi/dX and Gj^ = dGi/dcp. We mention that the above Lagrangian 
was first discovered by Horndeski [5] in a different but equivalent form [32]. 

Apart from the above scalar-tensor sectors, in a realistic cosmological scenario one 
needs to take into account the matter content of the universe, described by the Lagrangian 
C m , corresponding a perfect fluid with energy density p m and pressure p m . Then the total 
action is given by 

S = [ d 4 xV=g (C + C m ) , (2.6) 



where g is the determinant of the metric g^. 

In the following we impose a flat Friedmann-Robertson- Walker (FRW) background 
metric of the form ds 2 = —N 2 (t)dt 2 + a 2 (t)dx 2 , where t is the cosmic time, x l are the co- 
moving spatial coordinates, N(t) is the lapse function, and a(t) is the scale factor. Varying 
the action (2.6) with respect to N(t) and a(t) respectively, and setting N = 1, we obtain 

2XK iX -K + 6X<j)HG 3>x - 2XG 3 ,<p - 6H 2 G 4 + 2AH 2 X{G^ X + XG 4 ,xx) 
-l2HX<j>G 4><f>x - 6i^G 4 , + 2H 3 X<P(5G 5 , X + 2XG 5 , X x) 

-6H 2 X (3G 5 , + 2XG, )y4 ,x) = - Pm , (2.7) 
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K - 2X(G 3 ^ + 4>G 3) x) + 2(3H 2 + 2H)G 4 - 12H 2 XG^ X - 4HXG±,x - 8HXG ijX 
-8HXXG 4>XX + 2(4> + 2H<j))G^ + AXG^ H + 4X(4> - 2H4>)G^ X 
-2X{2H'^ + 2HH<j> + 3H 2 4>)G 5 , X - 4H 2 X 2 4> G 5t xx + <±HX(X - HX)G^ 4>X 
+2[2{HX + HX) + 3H 2 X]G 5)4> + 4HX0 G bM = -p m , (2.8) 

where dots denote derivatives with respect to t, and we also defined the Hubble parameter 
H = a I a. Variation of (2.6) with respect to (p(t) provides its evolution equation 



1 d 
a?~dl 



a 



j)=P (j> , (2.9) 



with 



J = <j)K )X + 6HXG 3)X - 2<PG 3 ^ + 6H 2 <j>{G 4)X + 2XG 4 y xx) - l2HXG^x 

+2H 3 X(3G 5 , X + 2XG 5:X x) - 6H 2 <p(G 5 ^ + XG^x) , (2.10) 
= Kj - 2X (G 3M + 4>G 3 ^ x ) + 6(2H 2 + H)G^ + 6H(X + 2HX)G^ X 
-6H 2 XG 5M + 2H 3 X<j> G^x • (2.11) 

Finally, the evolution equation for matter takes the standard form 

p m + 3H( Pm +p m ) = 0. (2.12) 

We mention here that the four equations (2.7), (2.8), (2.9) and (2.12), are not independent 
due to the Bianchi identities. In particular, the scalar equation (2.9) can be acquired from 
the other three equations [18]. 

In order to be able to perform the dynamical analysis of Galileon cosmology we need 
to focus on some more specific models. One class of Galileon scenarios has the above 
Lagrangian with the ansatzes: 

K(<f>,X) = X- V{4>), G 3 (cf>,X) = -gitfX, G A (<p,X) = G 6 (<f>,X) = 0, (2.13) 



corresponding to the action 



/ 



S = d xij—g 



R ~ \d^d^ - V{<t>) - Igi^d^d^ + C 



16ttG 2 r ^ ^ J 2 



(2.14) 



Such an action is able to capture the basic, quite general, and more interesting Galileon 
terms (one could straightforwardly include ansatzes with higher powers of X, such as the 
covariant Galileon model [7], however for simplicity we remain to the above simple but 
non-trivial Galileon action). In this case, the gravitational field equations (2.7) and (2.8) 
become 

H 2 = (pDE + Pm) i (2-15) 



3 

H = -4irG(p DE +PDE + Pm+Pm), (2.16) 
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where we have defined the effective dark energy sector with energy density and pressure 
respectively: 

Pde = y (l-QgH<P + g,^ 2 ) +V{ct>) , (2.17) 

12 , 

Pde = y{ 1 + 2 ^ + 9,^ 2 )-V^). (2.18) 

The scalar field equation (2.9) becomes 

(f> + 3H<p + 2g^ 2 4> + l -g M ^ - 3gH<p 2 - 6gHj>4> - 9gH 2 2 + V t4> = , (2.19) 

and we can immediately see that using (2. 17), (2. 18) it can be rewritten to the standard 
form 

p DE + 3H( PDE + p DE ) = 0. (2.20) 
Furthermore, we can define the dark energy equation-of-state parameter as 

_PDE f (l + 2# + ^ 2 )-^) 

wde = = -TT7 : • (^1) 

PDE ^(i-egHt + g^j+V^) 

One can clearly see that in this scenario, according to the form of g{4>), wde can be 
quintessence-like, phantom-like, or experience the phantom divide crossing during the evo- 
lution, which is a great advantage of Galileon cosmology. 

Without loss of generality in the following we restrict the analysis to the dust matter 
case, that is we assume that w m = p m /p m = 0. In this case it is convenient to introduce two 
additional quantities with great physical significance, namely the "total" equation-of-state 
parameter 

Pde 8ttG 



Wtot — — „ rj 2 

PDE + Pm 3H Z 

and the deceleration parameter 



^ (l + 2# + ^ 2 ) -V(<f>) 



(2.22) 



Q = - 1 ~§2=l + l w tot- (2-23) 

We close this section by mentioning that in order for the above scenario to be free of 
ghosts and Laplacian instabilities, and thus cosmologically viable, two conditions must be 
satisfied [18, 21, 61]. In particular, using the ansatzes (2.13) and units where k = 8ttG = 1, 
these write as [18] 

c| ^ 6WlH ~ 3wl ~ 6 ? = 6pm > 0, (2.24) 



for the avoidance of Laplacian instabilities associated with the scalar field propagation 
speed, and 



{4w 2 + 9wf ^ 



Q S = K ■ ' 1J >0, (2.25) 
3wi 
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for the absence of ghosts, where in our case 



w 2 



+ 2H, 



+ 9,, 



6Hg<p 



9H 2 



(2.26) 
(2.27) 



Finally, we stress that according to (2.21) and (2. 24), (2. 25) the phantom phase can be 
free of instabilities and thus cosmologically viable, as it was already shown for Galileon 
cosmology [18]. 



3 Phase space analysis 

In this section we perform a detailed phase-space and stability analysis of generalized 
Galileon cosmology. As usual, we first transform the above dynamical system into its au- 
tonomous form X' = f(X) [106-110], where X is the column vector constituted by suitably 
chosen auxiliary variables, f(X) the corresponding column vector of the autonomous equa- 
tions, with primes corresponding to derivatives with respect to M = In a. Next we extract 
its critical points X c demanding X'|x=x c = 0, and in order to determine their stability 
properties we expand around X c as X = X c + U, with U the column vector of the per- 
turbations. Therefore, for each critical point we expand the perturbation equations up to 
first order as U' = Q • U, where the matrix Q contains the coefficients of the perturbation 
equations. Lastly, the eigenvalues of Q calculated at each critical point determine its type 
and stability. 

In the scenario at hand we introduce the auxiliary variables: 

K<j) 



V6H 



z = g w H^> 
1 



(3.1) 



Using these variables the Friedmann equation (2.15) becomes 

/ n a \ 2 , 2 . VQzg'((f>)x 3 p m 
(1 — 6z)x + y + 



<?(0) 



+ 



3H 2 



1. 



(3.2) 



Moreover, using (3.2) and (2.17) we can write the matter and dark energy density param- 
eters as: 



n, 



n 



DE 



3H 2 



K 2 PDE 



3H 2 



9m 



(l-6z)x* + yi + V ~ 6Z9 ' {(t>)x3 



g{4>) 



(3.3) 
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Note that in the limit g((f>) —> the above quantities are well-defined, and they coincide 
with the usual quintessence ones [106]. 

In order to proceed forward we need to consider a specific scalar-field potential V{<j>) 
and a specific coupling function g(4>) with the Galileon term. Concerning V(4>) the usual 
assumption in dynamical investigations in the literature is to assume an exponential po- 
tential of the form 

V{cj>) = V e x ^, (3.4) 

since exponential potentials are known to be significant in various cosmological mod- 
els [106-112] (equivalently, but more generally, we could use potentials satisfying Ay = 
— v ^ ~ const, which is the case for arbitrary but nearly flat potentials [113-115]). 

Concerning g(<p), and in order to remain general, we will consider two ansatzes, namely 
the exponential one 

g(<P) = g e x ^, (3.5) 

and the power-law one 

g{4>) = go<P n - (3.6) 

The corresponding analysis will be performed separately in the following two subsections. 

3.1 Scenario 1: Exponential potential and exponential coupling function 

In this subsection we consider the exponential potential (3.4) and the exponential coupling 
function (3.5). In this case, using the auxiliary variables (3.1), the equations (2.15), (2.16) 
and (2.19) can be transformed to the autonomous form 



18z 2 x 3 + 4V6zX g x 2 + (2 - 12z)s 



1 {l8z 2 X 2 g x G + (9V6z\ g - 27V6z 2 X g | 



+ [54z 2 + (-6A 2 - 45) z + 3] x 4 + 6V6zX g - 3V6y 2 z(2X g + Ay) 
+ [(27 z - 3)y 2 - 9z - 3] x 2 - V6y 2 X v x} , 



(3.7) 



I8z 2 x 3 + iV^zXgX 2 + (2 - 12z)x 



{l8yz 2 X 2 g x 5 + yx 4 



9v / 6(A v - 2X g )z 2 + 9\/6A £ 



+ yx 3 [54z 2 + 12(A 9 A y - 3)z + 3] + x [(18z - 3)y 3 + (3 - 18%] 

+x 2 {y [6V6z{Xg - X v ) + \/6Ay] - 3V6y 3 z(2X g + Ay)}} , (3.8) 



18z 2 x 3 + A^QzXgX 2 + (2 - I2z)x 1 |-18z 3 A^ 5 + (27VQz 3 X g - 9\/6z 2 A, 
+ x 3 [-54z 3 + 9 (2X 2 g + 3) z 2 - 3z] + x [(3z - 9z 2 ) y 2 + 27z 2 - 9z] 



+x' 



-mVQXgZ 2 + 3V6y 2 (2X g + X v )z 2 + 2V6X g z - V6y 2 zX v ^ , (3.9) 



defined in the (non-compact) phase space ^ = { (x, y, z) : y > 0, z (y6xX g — 6) x 2 + x 2 + y 2 < l} 
Note that in this case the variable v is not needed. 
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Using (2.17) and the Friedmann equation (3.2), we can write the density parameters 



as: 

P 



o = < " l - i 
m ~ 3iP " 

K 2 PDE 



(1 - Qz)x 2 + y 2 + V6zx 3 X 



(3.12) 



q de = -^Sr = C 1 - 6z ) x2 + y 2 + vWA 9 , (3.10) 

O-TZ 

while for the dark-energy equation-of-state parameter (2.21) we obtain: 

6z 2 A 2 x 4 + 3\/6(l - 2z)zA 5 x 3 + [3z(3z - 4) + l]a; 2 - y/&y 2 z{2X g + X v )x + y 2 (6z - 1) 

UJDP = 

[zi 2 (v/6xA fl - 6) + + y 2 ] [z (9zx 2 + 2\/6A 9 x - 6) + l] 

(3.11) 

As we mentioned above, according to the variable z, that is according to the coupling 
function g((j>), in this scenario wde can be quintessence-like, phantom-like, or experience 
the phantom divide crossing during the evolution. Furthermore, the total equation-of-state 
parameter (2.22) becomes 

- PDE n 
Wtot = ■ = WdeMde = 

Pde + Pm 

6z 2 X 2 g x 4 + 3^/6(1 - 2z)zX g x 3 + [3z(3z - 4) + l]a; 2 - y/6y 2 z(2X g + \ v )x + y 2 (6z - 1) 

z (9zx 2 + 2V6X g x - 6) + 1 

and the deceleration parameter (2.23) reads 

H 1 3 
q = -l- 1 ^ = - + -w tot 

= \lz (9zx 2 + 2V6X g x - e) + 2} |9\/6(1 - 2z)zX g x i + [36(z - l)z + 3]x 2 

+18z 2 X 2 g x A + V6z [(2 - 6y 2 ) X g - 3y 2 X v ] x + (3y 2 - l) (6z - 1)} . (3.13) 

Finally, the two instability-related quantities (2.24) and (2.25) are respectively written as 

4 = [x [z (9zx 2 + 2y/EX g x - 6) + 1] 2 } 1 {3V6z 3 X g x 4 + 3z 2 x 3 (2X 2 g - 6z + 5) 

-27z 4 x 5 + 2\/6(l - iz)zX g x 2 + {z [3z (5 - 3y 2 ) - 4] + l} x + V6y 2 zX v }(3.U) 

and 



3x 2 [z (9zx 2 + 2y/6X g x - 6) + l] 
(3zx 2 + 1) 



= LJ ~ „ . " 1 L - (3.15) 



3.1.1 Finite phase-space analysis 

Let us now proceed to the phase-space analysis. The real and physically meaningful crit- 
ical points (x c ,y c ,z c ) of the autonomous system (3.7)-(3.9) (that is corresponding to an 
expanding universe, and thus possessing H > 0, with < £Ide < 1), are obtained by 
setting the left hand sides of the equations to zero, and they are presented in Table 1. 
In the same table we provide their existence conditions. For each critical point of Table 
1 we calculate the 3x3 matrix Q of the linearized perturbation equations of the system 
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Cr P 




1 1 

yc 




r ,"V1 TOT 


Sit a ni li tv 

k J Lex Ul±l L y 


A + 


+ 1 


o 


o 


always 


unstable for Ay > — V^Ag > 
saddle point otherwise 


A~ 


— 1 








always 


unstable for Ay < \/6,\g < 
saddle point otherwise 












unstable for X g > \/6,Ay > — A g 


B+ 


A s + a/ A g 2 — 6 





3-oT(A„) 


A^>6 


saddle point otherwise 


V6 


9 












unstable for A g < — \/6,Ay < — A g 


B~ 







3-a+(A g ) 


A*>6 


saddle point otherwise 


VD 


9 












stable node for — -\/3 < Ay < 0,A g < —Ay 


C 


X v 
V6 







< \\ < 6 


stable node for < Ay < \/3,A g > —Ay 










saddle point otherwise 


Co 





i 





Ay = 


stable node 












stable node for — < Ay < — \/3,A g < —Ay 


D 


V6 
2X V 


^6 
2X V 





Ay >3 


stable node for V3 < Ay < y^",A g > — Ay 
stable spiral for Ay < — y^",A g < — Ay 
stable spiral for Ay > y^,A g > —Ay 













always 


saddle 



Table 1. Scenario 1: Exponential potential and exponential coupling function. The real and 
physically meaningful critical points of the autonomous system (3.7)-(3.9), and their existence and 
stability conditions. We have introduced the functions a ± (A 9 ) = A g ± X g ^J X 2 g — 6. 

(3.7)-(3.9), and in order to determine the type and stability of the point we examine the 
sign of the real part of the eigenvalues of Q. The details of the analysis and the various 
eigenvalues are presented in Appendix A.l, and in Table 1 we summarize the stability 
results. Moreover, for each critical point we calculate the values of Qde, wdb, wtot and 
q given by (3.10)-(3.13), as well as the instability-related quantities c| and Qs given in 
(3. 14), (3. 15), and we summarize the results in Table 2. 

3.1.2 Phase-space analysis at infinity 

Since the dynamical system (3.7)-(3.9) is non-compact, there could be features in the 
asymptotic regime which are non-trivial for the global dynamics. Therefore, in order to 
complete the phase-space analysis we have to extend the investigation using the Poincare 
central projection method [116]. 

We consider the Poincare variables 

x r = pcosO sintp, z r = ps'mO sinip, y r = pcosijj, (3.16) 

with p = -E=g , t = \J x 2 + y 2 + z 2 , 9 G [0, 2n], and — f < ^ < § (we restrict the tp angle 
to this range since the physical region is given by y > 0) [116-119]. Therefore, the points 
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Cr. P. 




WDE 


Wtot 


q 


r 2 


Qs 


A± 


1 


1 


1 


2 


1 


3 


B ± 


1 


1 + a± ™ 

O 


1 + Q >) 

o 


i + a± ^ 





arbitrary 


C 


1 


X 2 

-1 + 

3 


X 2 
3 


X 2 
1 2 


1 


2 


Co 


1 


-1 


-1 


-1 


1 





D 


3 

\ 2 
•V 








1 

2 


1 


9 
2A^ 


Oi 





arbitrary 





1 

2 


1 






Table 2. Scenario 1: Exponential potential and exponential coupling function. The real and 
physically meaningful critical points of the autonomous system (3.7)-(3.9), and the corresponding 
values of the dark-energy density parameter £Ide , of the dark-energy equation-of-state parameter 
wde, of the total equation-of-state parameter Wtot and of the deceleration parameter q. The last 
two columns contain the instability-related parameters c| and Qs (from (3. 14), (3. 15)), which must 
respectively be non-negative and positive for a scenario free of ghosts and instabilities. We have 
introduced the functions a ± (A 9 ) = ± \ g J\ 2 g — 6 > 0. 



at "infinity" (r — > +oo) are those with p — )■ 1. Moreover, the physical phase space is now 
given by (x r ,y r ,z r ) G [—1,1] x [0,1] x [—1,1], such that 

x 2 + y 2 T x 2 r z r (V6X g x r - 6x 2 - 6y 2 - 6z 2 + 6) < ^ 



\-x 2 -y 2 - z 2 (i _ x 2 _ y 2 _ z 2f2 

and x 2 + y 2 + z 2 < 1. 

Inverting relations (3.16) and substituting into (3. 10), (3. 11), we obtain the dark energy 
density and equation-of-state parameters as a function of the Poincare variables, namely: 

x 2 + y 2 x 2 z r (V6X g x r - 6x 2 - 6y 2 - 6z 2 + 6) 

^ DE ~ l- x 2_ v 2_ z 2 ~ 2 ~2 2^ ' ^> 

i x r y r z r [i _ x z _ yz _ z z^ i 



WDE 



with £ 



V6z r X g x 3 r + (( - 6z r ^/() x 2 r + y 2 ( (9x 2 z 2 - 6( 3/2 z r + 2V6x r (\ g z r + C 2 ) } 1 
C {6z 2 \ 2 g xt + sVGzr (VC ~ 2j8r) \/C\ g x 3 r + ( (Qz 2 - 12^/Qz r + C) x 2 

-V6y 2 .z r ((2\ g + \ v )x r + y 2 (dz r - V?) C 3/2 } , (3.19) 

2 , and similarly substituting into (3.12), (3.13) we obtain the 
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Cr. P. 




Vr 




Stability 


&-DE 


WDE 


Wtot 


q 


r 2 
C S 


Qs 










±1 


unstable 





arbitrary 





i 

2 


arbitrary 






Table 3. Scenario 1: Exponential potential and exponential coupling function. The real and 
physically meaningful critical points at infinity of the autonomous system (3.7)-(3.9), and the cor- 
responding values of the dark-energy density parameter Qde, of the dark-energy equation-of-state 
parameter wde, of the total equation-of-state parameter wtot and of the deceleration parameter 
q. The last two columns contain the instability-related parameters c| and Qs (from (3. 22), (3. 23)), 
which must respectively be non-negative and positive for a scenario free of ghosts and instabilities. 



corresponding expressions for the total equation-of-state and deceleration parameters: 



Wtot 



C [9x 2 ,z 2 - 6( 3/2 z r + 2V6x r (\ g z r + C 2 )J [3V6z r ^ ~ 2z r J y/(X g x 3 r 



+C (9z 2 - Yl^lz r + C) x 2 r - y/(iy?z r ((2\ g + X v ) 
+6z 2 X 2 xf + y 2 (&z r - \/C 



(3.20) 



2( (9x 2 z 2 - 6( 3/2 z r + 2V6x r (X g z r + ( 2 ) J ~ j9\/fe r (y^ ~ 2z r J y^XgX 3 
+3C (l2z 2 - Ylyflzr + C) x 2 + V6z r ( [2(X g - 3y 2 (2X g + Ay)] x r 

+18z r 2 A^ + C 3/2 (yC-6z r ) (C-3y r 2 )}. (3.21) 

Finally, substitution into (3. 14), (3. 15) provides the corresponding expressions for the instability- 
related quantities: 



4 



x r (9x 2 z 2 . - 6( 3/2 z r + 2V6x r (X g z r + C 



-27zfx 5 r + 3V6z 3 (X g xf 



+3z 2 ( 3/2 {2X 2 g + 5) - 6z r ] x 3 + 2V6z r ( y/l - 4z r ) C 5/2 A 



+c 2 



C (l5z r 2 - Ay/CZr + C 



9y r z r 



? r + \/6y 2 z r C 3 Ay}, (3.22) 



3x 2 (9x 2 z 2 - 6( 3 / 2 z r + 2V6x r (X g z r + ( 2 ) 
{3z T x? + C 3 / 2 ) 2 



(3.23) 



Therefore, performing the analysis described in Appendix A. 2 we conclude that there 
are two physical critical points at infinity, namely K . These critical points, along with 
their stability conditions and the corresponding values of observables, are presented in 
Table 3. 

3.2 Scenario 2: Exponential potential and power-law coupling function 

Let us consider an exponential potential V(<p) = Vae^^ and a power-law coupling function 
= go<p n - In this case, using the auxiliary variables (3.1), the equations (2.15), (2.16) 
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and (2.19) can be transformed to the autonomous form: 

x' = 2V6 + 2Axvzn-12zV6 + l8V6x 2 z 2 * |l8\/6n (n + 1) v 2 z 2 x 5 - 6y 2 \ v 
+Mvznx 4 (1 - 3z) + 3\/6x 3 {l8z 2 + [2n(l - n)v 2 - 15] z + l} 
-18zx 2 [(2rm + Ay) y 2 - 2vn] + 3V&c [(9z - 1) y 2 - 3z - l] } (3.24) 



2^6 + 24x?;zn - 12z\/6 + ISV&eV 1 18\/6n (n + 1) v 2 z 2 yx 4 

+3\/6yx 2 [l8z 2 + 4 (Aym; - 3) z + l] + 3\/6 (6z - 1) y 3 
— 6x {3zy 3 (2m; + Ay) — y [6z (nu — Ay) + Ay]} 
+3v / 6y (1 - 6z) + 54yx 3 [(Ay - 2nv) z 2 + vzn] } 

2x (y/E + \2xvzn - 6zVE + 9\/6x 2 z 2 ) 1 |-18\/6n (n + 1) v 2 z 3 x 5 

+54nvz 2 x A (3z - 1) - 3\/(3zx 3 {l8z 2 - [2n (1 + 3n) ^ 2 + 9] z + l} 
+6zx 2 [(6m + 3Ay) zy 2 — 18nvz + 2i>n] 

+3«>/6a! [y 2 (1 - 3z) + 9z - 3] - 6zy 2 Ay} 



t/ = — y/Ev 2 x, 
defined in the non-compact phase space 



(3.25) 



(3.26) 
(3.27) 



[x, y,z,v) : y > 0, z (y/Envx — 6^ x 2 + x 2 + y 2 < 1, t> € k| 



Note that contrary to the previous Scenario 1, in the case at hand we do need the fourth 
auxiliary variable v. 

Using (3.2) and (2.17) we can write the density parameters as: 

Pm 



n 



DE 



3H 2 ~ 
_ fc 2 p<f> 
" 3# 2 



1 



(1 - 6z)x 2 + y 2 + y/& 



nvzx 



(1 - 6z)x 2 + y 2 + y/6: 



nvzx 



(3.28) 



while for the dark-energy equation-of-state parameter (2.21) we obtain: 

wde = | (Vftx 2 — gV6zx 2 + 6vx 3 zn + v^y 2 ^ + \2xvzn — 6zV6 + 9v / 6x 2 z 2 ^ j 
6 6 (1 + n) nv 2 z 2 x 4 + 3\/6n (1 - 2z) vx 3 z + (l + 9z 2 - 12 z) x 2 

-V6 (2nv + 6X V ) zy 2 x + (6z - 1) y 2 . (3.29) 
Furthermore, the total equation-of-state parameter (2.22) parameter reads 



Wtot 



1 + 2y/h~xvzn — 6z + 9x 2 z 2 \ 6 (1 + n) nv 2 z 2 x 4 + 3\/6n (1 — 2z) 



vx z 



+ (1 + 9z 2 - 12z) x 2 - V6 {2nv + 6Ay) zy 2 x + (6z - 1) y 2 j , (3.30) 
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and the deceleration parameter (2.23) becomes 



1 3 

2 + 2 WM - 



(3.31) 



Finally, from (2. 24), (2. 25) we find 

1 + 2VE zxnv 



4 



6z + 9x 2 z 2 



-27x 5 z 4 + 



nv 



+3z 2 x 3 [-6z + 5 + 2n (n - 1) v 2 ] + 2y/6znv (1 - Az) x 2 
+x (l5z 2 - 9y 2 z 2 - Az + l) + z\/oy 2 Ay} 



and 



(3 + 6V6zxnv - lBz + 27xV) 
(3x 2 z + l) 2 



(3.32) 



(3.33) 



3.2.1 Finite phase-space analysis 

Let us now proceed to the phase-space analysis. The real and physically meaningful critical 
points (x c , y c , z c , v c ) of the autonomous system (3.24)-(3.27) (that is corresponding to an 
expanding universe, and thus possessing H > 0, with < floE < 1)> are presented 
in Table 4, along with their existence conditions. The details of the analysis and the 
various eigenvalues are presented in Appendix B.l. Furthermore, in Table 5 we display the 
corresponding values of the observables Qde, wde, wtot an d q given by (3.28)-(3.31), as 
well as the instability-related quantities Cg and Qs given in (3. 32), (3. 33). We mention 
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always 
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stable (not asymptotically for n^O) 
stable (asymptotically for n = 0) 














stable node for 3 < Ay < ^ 
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2\ v 


2\ v 








X 2 y>3 


stable spiral for < 
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always 


saddle 



Table 4. Scenario 2: Exponential potential and power-law coupling function. The real and phys- 
ically meaningful critical points of the autonomous system (3.24)- (3.27), and their existence and 
stability conditions. 

that the stability of the above points does not depend on the exponent n, since these points 
have z c = in which case n disappears from the equations and v decouples. Similarly, since 
they correspond to v = 0, that is to <j) — > 00, their coordinates themselves do not depend 
on n (in other words for (ft — > 00 all exponents n 7^ 0, of the same sign, are equivalent). 
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Table 5. Scenario 2: Exponential potential and power-law coupling function. The real and phys- 
ically meaningful critical points of the autonomous system (3.24)- (3.27), and the corresponding 
values of the dark-energy density parameter £Ide, of the dark-energy equation-of-state parameter 
Wde, of the total equation-of-state parameter Wtot and of the deceleration parameter q. The last 
two columns contain the instability-related parameters c s and Qs (from (3. 32), (3. 33)), which must 
respectively be non-negative and positive for a scenario free of ghosts and instabilities. 



3.2.2 Phase-space analysis at infinity 

Due to the fact that the dynamical system (3.24)-(3.27) is non-compact, there could be 
features in the asymptotic regime which are non-trivial for the global dynamics. Thus, 
in order to complete the analysis of the phase space we must extend our study using the 
Poincare central projection method [116]. 
We consider the Poincare variables 

x r = p cos 9 sin <3? sin ip,y r = p cos tp, z r = p sin 9 sin $ sin tp, v r = p cos $ sin ip, (3.34) 



where p = —^==^, r = \J x 2 + y 2 + z 2 + v 2 , #,<3? G [0, 27r], and — | < ip < ^ (we restrict 
the angle ip to this interval since the physical region is given by y > 0) [116-119]. Thus, 
the points at infinity (r — > +00) are those having p — )■ 1. Furthermore, the physical phase 
space is given by (x r ,y r , z r ,v r ) G [—1,1] x [0,1] x [—1,1] x [—1,1], such that 

x 2 r + yl x 2 z r (\/6nv r x r - 6v 2 - Gxj - Gy 2 - 6z 2 + 6) < ^ ^ 



v 2 -x 2 -y 2 - z 2 (x_ v 2_ x 2_ y 2 _ z 2f/z 



and v 2 + x 2 + y 2 . + z 2 < 1. 

Inverting relations (3.34) and substituting into (3. 28), (3. 29), we obtain the dark energy 
density and equation-of-state parameters as a function of the Poincare variables, namely: 

x 2 + y 2 x 2 z r (^/6X g x r - 6x 2 - 6y 2 . - 6z 2 + 6) 

= Z o o o (o.OD) 

l-x 2 -y 2 .- z 2 (l-x 2 - y 2 - z 2 f 2 
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WDE 



2z r y/l (V6nv r x r - 3C) + 9x 2 z r 2 + ( 2 ] [V6i 



; r X^Z r + C 3 ^ 2 (^r + Vr) ~ Gx^ZrC 



r 



C 3 / 2 (3x r 2 - 2 y 2 ) 



a/C 6n(n + l)f 2 x^2 2 — Gv^mv^r^rC + V6nv r x r z r 



+6z r ( 5/2 (y 2 - 2x1) + XrZrC (Vxrzr - VGy 2 \) + C 3 ( 



(3.37) 



where C = 1 



2 — y 2 . — z 2 . — v 2 , and similarly substituting into (3.30), (3.31) we obtain 
the corresponding expressions for the total equation-of-state and deceleration parameters: 



W t ot 



C 2 (2V6y^nv r x r z r + 9x 2 z 2 - 6z r ( 3/2 + C 2 )l [6n(n + l)v 2 .xjz 2 



nv r x^ r z 2 



( + v / 6nu r a; r .^C 3/2 (3x 2 - 2y 2 ) + 6z r ( 5 / 2 (y 2 - 2x 2 ) 

+x r z r ( 2 (9x r z r - V6y 2 \) + ( 3 (x 2 - y 2 )] , (3.38) 



q = {2C 2 2z rA /C (VfiravaV - 3C) + 9s 2 2 2 + C 2 } {l8n(n + V)v 2 r x^ r z 2 r 

+9V6nv r xfz r ( (y/c - 2z r ) + V6x r z r C 3/2 2nu r (C - 3y 2 ) - 3y 2 \/CA 

+3x 2 C 2 (l2z 2 - I2z r y^+C) + C 5/2 (C " 3y r 2 ) (^C - 6z r ) } .(3.39) 

Finally, substitution into (3. 32), (3. 33) provides the corresponding expressions for the instability- 
related quantities: 

-1 



4 



3C ) + 9x 2 z 2 + C 



3x r z r £ 



2(n — l)m; 2 — 6z r \/C + 5C + Sv^mva^rv/C — 27x: 



2\J~Qnv T x 2 z r C, 2 ( -y/C — 4z, 

4_3 



5 r zj 



C (l5z 2 - 4z r + C) - 9y r 2 z 2 ] + v / 6 2 / 2 z r C 3 A} , (3.40) 



Qt 



3x 2 [2z r y/( (V6nv r x r - 3C) + 9x 2 z r 2 + ( 2 
(3x 2 z r + C 3 / 2 ) 2 



(3.41) 
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Table 6. Scenario 2: Exponential potential and power-law coupling function. The real and phys- 
ically meaningful critical points at infinity of the autonomous system (3.24)-(3.27), and the corre- 
sponding values of the dark-energy density parameter Qde, of the dark-energy equation-of-state 
parameter wde, of the total equation-of-state parameter wtot and of the deceleration parameter 
q. The last two columns contain the instability-related parameters c| and Qs (from (3. 40), (3. 41)), 
which must respectively be non-negative and positive for a scenario free of ghosts and instabilities. 

Therefore, performing the analysis described in Appendix B.2 we are led to the critical 
points at infinity and M^, which are presented in Table 6 along with their stability 
conditions. 
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4 Cosmological Implications 



Since we have performed the complete phase-space analysis of a quite general subclass of 
generalized Galileon scenario, we can now discuss the corresponding cosmological behavior. 
In the following subsections we analyze the two different scenarios, namely that of expo- 
nential potential and exponential coupling function, and that of exponential potential and 
power-law coupling function, separately. In each scenario we analyze the observables at the 
solutions, and moreover the instability-related quantities cs and Qs, examining whether 
they satisfy (2. 24), (2. 25) in order for the corresponding solution to be physical. However, 
we mention that this ghost investigation is applied at the critical points, thus even if they 
are found to be ghost and instabilities free, it is not a proof that the universe evolution to- 
wards them did not pass through a state with ghosts and Laplacian instabilities. Therefore 
in order to obtain a physical evolution of the universe one has to ensure that conditions 
(2. 24), (2. 25) are satisfied everywhere in the examined evolution, and not only at late times. 

4.1 Scenario 1: Exponential potential and exponential coupling function 

The critical points of the scenario at hand are the following. Points exist always, that 
is for every values of the scenario parameters Ay,A 9 ,Vo,go> they are unstable or saddle, 
and thus they cannot be the late-time state of the universe. They correspond to a non- 
accelerating, dark-energy dominated universe, with a stiff dark-energy equation-of-state 
parameter equal to 1. Finally, the instability-related quantities cs and Qs do satisfy the 
corresponding conditions (2. 24), (2. 25), namely cs > and Qs > 0, and thus these solutions 
are free of instabilities. Both of them exist in standard quintessence [106]. 

Point 0\ is a saddle one and thus it cannot attract the universe at late times. It cor- 
responds to a non-accelerating, dark-matter dominated universe, with zero total equation- 
of-state parameter. The instability-related quantities cs and Qs satisfy (2. 24), (2. 25), and 
thus this solution is free of instabilities (the fact that is exactly and not positive is 
not a problem, since this happens only at one point and it is not zero identically, in which 
case there might be a non-perturbative ghost [120, 121] 1 ). This point exists in standard 
quintessence too [106]. 

Point C exists for < \y < 6 and it is a stable one in the region of the parameter 
space shown in Table 1, and thus it can attract the universe at late times. It corresponds 
to a dark-energy dominated universe, with a dark-energy equation-of-state parameter lying 
in the quintessence regime, which can be accelerating or not according to the Ay-value. 
Additionally, this solution is free of instabilities. This point exists in standard quintessence 
[106]. It is quite important, since it is both stable and possesses wde and q compatible 
with observations. 

Furthermore, we mention that in the specific case where Ay = 0, that is in the case of 
constant or zero usual potential, there exist the critical point Co, and it is always stable. 
It corresponds to the de Sitter solution, where the universe is accelerating and dark-energy 
dominated, with the dark-energy behaving like a cosmological constant (wde = — 1) 5 and 
it is free of instabilities. This point exists in standard quintessence [106]. It is quite 

1 We thank A. De Felice for this comment. 
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important, since it is both stable and possesses wde and q compatible with observations. 
Although at first sight it seems to be the Ay — > limit of point C, this is not the case 
since the eigenvalues of Co do not arise from the Ay — > limit of C-eigenvalues (Ay = is 
a bifurcation value) and thus it has to be considered separately. Moreover, since in many 
Galileon works the authors do not consider a potential, point Co just gives the late-time 
state of the universe in these cases [18, 61]. 

Point D exists for X v > 3 and in this case it is always stable, that is it can be the 
late-time state of the universe, and it is free of instabilities. It has the advantage that the 
dark-energy density parameter is in the interval < £Ide < 1 5 that is it can alleviate the 
coincidence problem since dark energy and dark matter density parameters can be of the 
same order. However, it has the disadvantage that it is not accelerating and wde = 0, 
which are not favored by observations. This point exists in standard quintessence [106] 
too. 




xr 



Figure 1. Trajectories in the y r -x r plane of the Poincare phase space for the Scenario 1, that 
is for exponential potential and exponential coupling function. We use Ay = 1 and X g arbitrary 
(for the numerics we choose X g = 1 but different X g 's correspond to the same projection on y r -x r 
plane). The region inside the inner semi-circle (seen as semi-ellipse in the figure scale), marked by 
the thick dashed-dotted line, is the physical part of the phase space. The region above the dotted line 
marks the region corresponding to accelerating universe (q < 0). In this projection the dark-energy 
dominated, accelerating, quintessence-like solution C is a stable solution, 0\ is saddle point, and 
A are unstable. 



Apart from the above points that exist also in standard quintessence, the scenario at 
hand possesses two additional critical points, namely B^. They correspond to dark-energy 
domination, with a dark-energy equation-of-state parameter lying in the quintessence 
regime, where the universe is non-accelerating (q > 0), and they are free of instabili- 
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ties. However, these points are not stable and thus they cannot attract the universe at late 
times. 

Finally, the present scenario possesses two critical point at infinity, namely K^. They 
correspond to a dark- matter dominated, non-accelerating universe, with arbitrary wde but 
with a zero total equation-of-state parameter wtot, which are also free of instabilities. They 
are always unstable and therefore they cannot be the late-time state of the universe. 

From the above analysis we observe that at the stable critical points, C and D, we have 
<fi — > 0, (j) — > — sign(A\/)oo and thus for XyX g > we obtain g(<f>) — > while for XyX g < 
we obtain g(4>) —> oo (for Ay = g{4>) can be zero, finite, or infinity). Similarly, for Co we 
see that for X g > we obtain g(<f>) — > while for X g < we obtain <?(</>) —> oo. In all cases, 
if X g = then obviously g(<j)) = const. From these we deduce that Galileons (simple or 
generalized) may survive at late-time cosmology or may be completely disappeared by the 
dynamics, depending on the model parameters. However, firstly we observe that Galileon 
cosmology possesses the same stable late-time solutions with standard quintessence [106]. 
Secondly, and more interestingly, as we observe from Table 2, the corresponding observables 
of these solutions do not depend on the Galileon terms, but only on the usual terms and 
especially on the standard scalar potential (note that even the instability-related quantities 
do not depend on the Galileon terms either). This is a main result of the present work. 




Figure 2. Trajectories in the y r -x r plane of the Poincare phase space for the Scenario 1, that is 
for exponential potential and exponential coupling function. We use \y = 2 and X g arbitrary (for 
the numerics we choose X g = 1 but different X g 's correspond to the same projection on y r -x r plane). 
The region inside the inner semi-circle (seen as semi-ellipse in the figure scale), marked by the thick 
dashed-dotted line, is the physical part of the phase space. The region above the dotted line marks 
the region corresponding to accelerating universe (q < 0). In this projection the non- accelerating, 
dust-like (wde = 0) solution D is a stable spiral, C and 0\ are saddle points, and A^ 1 are unstable. 
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Therefore, although the Galileon terms of generalized Galileon cosmology can play 
an important role at early (inflationary) times [28-35] or at recent times [18-27], in the 
future they will not have any effect on the universe evolution, when it will asymptotically 
reach its stable state. This is consistent with observational constraints, which disfavor the 
Galileon presence in the current universe [43, 47, 65-67, 89] (see also [70]). Finally, one 
could still ask whether going beyond the background analysis, and examining perturbation- 
related observables at the late-time stable solutions, would reveal a dependence of these 
observables on the Galileons. Although this complicated project lies beyond the scope of 
the present work, the examination of the instability-related quantities cs and Qs, which 
arise from a perturbative analysis, and their independence on the Galileon terms, is an 
indication that even the perturbation observables will not depend on the Galileon terms. 




Figure 3. Trajectories in the y r -x r plane of the Poincare phase space for the Scenario 1, that 
is for exponential potential and exponential coupling function, for the specific case Ay = and X g 
arbitrary (for the numerics we choose X g = 1 but different X g 's correspond to the same projection 
on y r -x r plane). The region inside the inner semi- circle (seen as semi- ellipse in the figure scale), 
marked by the thick dashed- dotted line, is the physical part of the phase space. The region above the 
dotted line marks the region corresponding to accelerating universe (q < 0). In this projection the 
de Sitter solution Cq is a stable node, 0\ is saddle point, and are unstable. 



In order to present the above behavior in a more transparent way, we first evolve the 
autonomous system (3.7)-(3.9) numerically for the choice Ay = 1 and X g arbitrary (with 
X g > — Xy), and in Fig. 1 we depict the corresponding phase-space behavior, projected in 
the Poincare x r -y r plane (for the numerics we choose X g = 1 but different A 9 's correspond to 
the same projection on y T -x r plane). As we can see, in this case the dark-energy dominated, 
accelerating, quintessence-like critical point C is the late-time solution of the universe. 

Similarly, in Fig. 2 we depict the corresponding phase-space behavior, for Ay = 2 and 
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Xg arbitrary (we use X g = 1, which satisfies X g > —Ay). In this case the non-accelerating 
critical point D is the late-time solution of the universe, while C is saddle point. Finally, 
in Fig. 3 we depict the phase-space behavior for the specific case Ay = and X g arbitrary. 
In this case the universe is attracted by the de Sitter solution Co as expected. 

4.2 Scenario 2: Exponential potential and power-law coupling function 

This scenario possesses the following critical points. exist for every value of the param- 
eters Xv,n,Vo,go, they are saddle points, and thus they cannot be the late-time state of the 
universe. They correspond to a non-accelerating, dark-energy dominated universe, with a 
stiff dark-energy equation-of-state parameter equal to 1, and they are free of instabilities. 
Both of them exist in standard quintessence [106]. 

Point O2 is a saddle one and thus it cannot attract the universe at late times. It cor- 
responds to a non-accelerating, dark-matter dominated universe, with zero total equation- 
of-state parameter, and it is free of instabilities. This point exists in standard quintessence 
too [106]. 

Point F exists for < Xy < 6 and it is a stable one for X v < 3 and thus it can be the 
late-time state of the universe. It corresponds to a dark-energy dominated universe, with 
a dark-energy equation-of-state parameter in the quintessence regime, which can be accel- 
erating or not according to the Ay-value. Additionally, this solution is free of instabilities. 
This point is quite important, since it is both stable and possesses wde and q compatible 
with observations. It exists in standard quintessence [106] too. 

There is another interesting critical point, namely Fq, which is obtained only in the case 
where Ay = 0, that is in the case of constant or zero usual potential, and it is always stable. 
It corresponds to the de Sitter solution, where the universe is accelerating and dark-energy 
dominated, with the dark-energy behaving like a cosmological constant {wde = — 1) ; and 
it is free of instabilities. At first sight it seems to be the Ay — > limit of point F, however 
this is not the case since the eigenvalues of -Fo do not arise from the Ay — > limit of 
-F-eigenvalues (Ay = is a bifurcation value) and thus it has to be considered separately. 
Furthermore, since in many Galileon works the standard potential is not considered, point 
Fq is straightforwardly the corresponding late-time state of the universe in these cases. 

Point G exists for Xy > 3 and in this case it is always stable, that is it can attract 
the universe at late times, and it is free of instabilities. It has the advantage that the 
dark-energy density parameter lies in the interval < £Ide < 1 5 that is it can alleviate the 
coincidence problem, but it has the disadvantage that it is not accelerating and possesses 
wde = 0, which are not favored by observations. This point exists in standard quintessence 
[106] too. 

Finally, the scenario at hand possesses the critical point at infinity and . 
They correspond to a dark-matter dominated, non-accelerating universe, with a zero total 
equation-of-state parameter wtoti which are also free of instabilities. are saddle points 
whereas are unstable, and thus they cannot be the late-time state of the universe. 

From the above analysis we can see that at the stable critical points, F and G, we have 
<\> — > 0, (j) — > — sign(Ay)oo and thus for n < we obtain g((p) — > while for n > we obtain 
g(4>) — > 00. Similarly, for Fq we see that for n < we obtain g((p) — > while for n > 
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we obtain g{4>) — > oo. In all cases, if n = then obviously g(4>) = const. Prom these it is 
implied that Galileons (simple or generalized) may survive at late-time cosmology or may 
be completely disappeared by the dynamics, depending on the model parameters. Similarly 
to the previous subsection, firstly we observe that this scenario possesses the same stable 
late-time solutions with standard quintessence [106]. Furthermore, as we observe from 
Table 5, the corresponding observables of these solutions do not depend on the Galileon 
terms, but only on the usual terms and especially on the standard scalar potential. Even 
the instability-related quantities do not depend on the Galileon terms either. This is a 
main result of the present work. 




Figure 4. Projection of orbits on the x r -y r -z r space of the Poincare phase space for the Scenario 
2, that is for exponential potential and power-law coupling function, for the specific case Ay = and 
n arbitrary (for the numerics we choose n = 1 but different n's correspond to the same projection). 
In this projection the de Sitter solution Fo is the attractor, whereas L~ and are unstable and 
O2 is a saddle point. 

Thus, and similarly to the scenario of the previous subsection, although the Galileons 
can play an important role at inflationary [28-35] or at recent times [18-27], in the fu- 
ture, when the universe will asymptotically reach its stable state, they will not have any 
effect on its evolution (this is consistent with observational constraints which disfavor the 
Galileon presence in the current universe [43, 47, 65-67, 89]. One could still ask whether 
the observables related to perturbations would depend on the Galileon terms, however the 
investigation of the instability-related quantities cs and Qs, which come from a perturba- 
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tive analysis, and their independence on the Galileon terms, is an indication that even the 
perturbation observables will not depend on the Galileon terms. 

In order to present the above behavior in a more transparent way, we first evolve the 
autonomous system (3.24)-(3.27) numerically for the choice Ay = and n arbitrary, and 
in Fig. 4 we depict the corresponding phase-space behavior, projected on the Poincare 
Xf-y^-Zf space (for the numerics we choose n = 1 but different re's correspond to the same 
projection). As we can wee, in this case the de Sitter critical point Fq is the late-time 
solution of the universe while L~ and are unstable. Finally, note that the sign of z r is 
invariant, thus the plane z r = cannot be crossed by the orbits. 

We close this section by comparing the present scenario of exponential potential and 
power-law coupling function, with the one of the previous subsection, that is with expo- 
nential potential and exponential coupling function. A first observation is that for X g = 
and n = the critical points of both scenarios coincide, which was expected since in this 
case the two models coincide and they both transit to the simple Galileon scenario, where 
the couplings are constants and not functions of the scalar field. Furthermore, in both sce- 
narios we found that there are not additional stable late-time solutions apart from those of 
standard quintessence. What is more interesting is that all the observables do not depend 
on the Galileon terms, that is the Galileons do not play any role at late times. However, 
in general the two scenarios are different, with the main difference being that they have 
different non-stable critical points and thus phase-space dynamics. 

5 Conclusions 

In the present work we have investigated the dynamical behavior of generalized Galileon 
cosmology, which is a recent construction based on higher derivatives in the action, along 
with the imposition of the Galilean symmetry, with the requirement of maintaining the 
equations of motion second ordered. However, one additionally extends the constant co- 
efficients of the various action-terms of the simple Galileon formulation, into arbitrary 
functions of the scalar field. 

Performing a detailed phase-space analysis of two distinct and quite general scenarios, 
namely of exponential usual potential and exponential Galileon coupling-function, and 
exponential usual potential and power-law Galileon coupling-function, we extracted the 
stable solutions, that is the solutions that will be the late-time states of the universe. In 
each of these stable late-time solutions we have calculated various observables, such as 
the dark-energy equation-of-state parameter, the deceleration parameter, and the dark- 
energy and matter density parameters. Furthermore, in order to examine whether these 
solutions are free of ghosts and instabilities, we calculate the corresponding instability- 
related quantities of the literature. 

In the case where the Galileon terms are absent we recovered the results of standard 
quintessence. In the case where the Galileons are present we found that at late times they 
may survive at or they may completely disappear by the dynamics, depending on the model 
parameters. However, independently of their disappearance or survival, we found that 
the scenario at hand possesses exactly the same stable late-time solutions with standard 
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quintessence, which are moreover free of ghosts and instabilities. More interestingly, the 
corresponding observables at these stable late-time solutions do not depend on the Galileon 
terms, but only on the usual action-terms and especially on the standard scalar potential. 
This is a main result of the present work, and it shows that although the Galileons can play 
an important role at inflationary [28-35] or at recent times [18-27], in the future, when 
the universe will asymptotically reach its stable state, they will not have any effect on its 
evolution. This is in agreement with observational constraints which disfavor the Galileon 
presence in the current universe [43, 47, 65-67, 89]. 

One could ask whether the above behavior is a result of the specific ansantzes, or it 
has a general character. Although one cannot exclude the case where suitably chosen or 
tuned ansantzes can lead to significant Galileon effects on the observables at late times, in 
general one expects the above behavior to be valid at late times in the large majority of 
cases, due to the downgrading of the role of the derivative (Galileon) terms in an eternally 
expanding universe. 

We close this work by mentioning that the obtained result that the Galileons will not 
play any role and that they will not have any observable effect at late times, was extracted 
for the background evolution and the corresponding observables. There could still be 
the case that Galileons could leave their signatures in observables related to perturbations. 
Although such an investigation is quite complicated and lies beyond the scope of the present 
work, the examination of instability-related quantities arising from perturbation analysis 
that we did perform in this work showed that they do not depend on the Galileon terms 
either, which is a quite strong indication that even the perturbation observables will not 
depend on the Galileon terms. 
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A Stability of Scenario 1 

A.l Stability of the finite critical points 

For the critical points (x c ,y c ,z c ) of the autonomous system system (3.7)-(3.9), the coeffi- 
cients of the perturbation equations form a 3 x 3 matrix Q, which can be easily obtained, 
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Table 7. The eigenvalues of the matrix Q of the perturbation equations of the autonomous system 
(3.7)-(3.9). We have defined the functions a ± (A ff ) = X g ± X g J\ 2 g - 6. 

however, since they are rather lengthy expressions we do not present them explicitly. De- 
spite this complicated form, using the specific critical points presented in Table 1 the 
matrix Q acquires a simple form that allows for an easy calculation of its eigenvalues. The 
corresponding eigenvalues Vi (i = 1,2,3) for each critical point are presented in Table 7. 

Thus, by determining the sign of the real parts of these eigenvalues, we can classify 
the corresponding critical point. In particular, if all the eigenvalues of a critical point have 
negative real parts then the corresponding point is stable, if they all have positive real 
parts then it is unstable, and if they change sign then it is a saddle point. 

We mention here that in order to examine correctly the stability of point Cq (since the 
second eigenvalue vanishes at first place) we introduce the local coordinates 

{x, y - 1, z] = e {x, y, z] + 0(e 2 ) 

where e is a constant satisfying eCl. Then we obtain the local evolution equations 

x' = — 35 + h.o.t, y = — 3y + h.o.t, z = — 3z + h.o.t, (A.l) 
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with h.o.t denoting "higher order terms". Prom the above linear approximation we deduce 
that the eigenvalues of the linearization for point Co are all negative and are displayed in 
Table 7, thus Co is a stable node. 

Similarly, in order to examine correctly the stability of point 0\ we introduce the local 
coordinates 

{x, y,z} = e {x, y, z] + 0(e 2 ) 

with e < 1. The local evolution equations valid around the origin then write as 

3 3 9 

x=—-x + h.o.t, y = -y + h.o.t, z! = — -z + h.o.t. (A. 2) 

The eigenvalues of the linearization for point 0\ are displayed in Table 7, and thus 0\ is 
a saddle point. 

A. 2 Stability of the critical points at infinity 

We consider the Poincare variables 

x r = pcos9 sinip, z r = ps'm9 s'mip, y r = pcosip, (A. 3) 



with p = -p==z, r = \J x 2 + y 2 + z 2 , 9 S [0, 2tt], and — \ < ip < \ (we restrict the ip angle 
to this range since the physical region is given by y > 0) [116-119]. Therefore, the points 
at infinity (r — > +oo) are those with p — > 1. Moreover, the physical phase space is now 
given by (x r ,y r ,z r ) G [—1, 1] x [0, 1] x [—1, 1] such that 

xl + y 2 x 2 z r (VGX g x r - 6xj - Gy 2 - Gz 2 + 6) 



1 - x 2 - y 2 - z 2 (i_ x 2_ y 2_ z 2f/i 

and x 2 + y 2 + z 2 . < 1. 

Performing the transformation (A. 3) on the system (3.7)-(3.9), and taking the limit 
p — > 1, the leading terms in the resulting system are 

p -> -X 2 g cos 2 6 sin 2 ip [cos(26») sin 2 ip + cos 2 ip] , (A.5) 
(1 - p 2 )6' -2A 2 sin 6 cos 3 6 sin 2 V, (A.6) 
(1 - p 2 )i>' -2A 2 sin 2 9 cos 2 9 sin 3 V cos tp. (A.7) 

Since the equation for p decouples, we only need to investigate the subsystem of the angular 
variables. Thus, equating (A.6) and (A.7) to zero, we obtain three classes of solutions: 
Class A: 

x T = sin-i/', y r = cost/', z r = 0, (A. 8) 

for some tp. These points lead to a divergence in £Ide in (3.18) and thus we deduce that 
they are unphysical. 
Class B: 

x r = 0, y r = cos?/', z r = ±sin^, (A.9) 
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for some ip. These points correspond to 6 — -^Y substituting these values for x r ,y r ,z r 
in the evolution equations for x' r , y' r , z' r we obtain that y r must be zero at a fixed point at 
infinity. Thus, we have the points at infinity given by : [x r = 0, y r = 0, z r = ±1] . 
Class C: 

x r = 0, y r = 1, z r = 0. (A. 10) 

This point leads to a divergence in SIde m (3.18) and thus it is unphysical. 

In summary, the scenario of exponential potential and exponential coupling function 
has only two critical points at infinity, namely : [x r = 0, y r = 0,z r = ±1]. In order to 
calculate the eigenvalues of the linearization of the dynamical system for x r ,y r ,z r around 

, we proceed as follows: having calculated the linear perturbation matrix Q we evaluate 
it at y r = 0, then we take the limit x r — > and finally we incorporate the lateral limit 
z r — > 1~ (that is from below) for K + or z r — > — 1 + (that is from above) for K~ . Calculating 
the eigenvalues of the resulting matrix we obtain that the eigenvalues of the linearization 
for both K± are given by (3,f ,§). Therefore, we conclude that K are unstable. 

B Stability of Scenario 2 

B.l Stability of the finite critical points 

For the critical points (x c , y c , z c , v c ) of the autonomous system (3.24)-(3.27), the coefficients 
of the perturbation equations form a 4 x 4 matrix Q. The explicit expressions for these 
coefficients are quite lengthy and for simplicity we do not present them explicitly. However, 
using the explicit critical points presented in Table 4 we can straightforwardly see that the 
matrix Q acquires a simple form that allows for an easy calculation of its eigenvalues. The 
corresponding eigenvalues (i = 1, 2, 3, 4) for each critical point are presented in Table 8. 
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Table 8. The eigenvalues of the matrix Q of the perturbation equations of the autonomous system 
(3.24)-(3.27). 
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Note that the stability of the above points does not depend on the exponent n, since 
these points have z c = in which case n disappears from the equations and v decouples. 
Similarly, since they correspond to v = 0, that is to <p ~ > 00 > their coordinates themselves 
do not depend on n (in other words for <f> —¥ 00 all exponents n ^ of the same sign are 
equivalent). 

Strictly speaking, since the eigenvalue of the fourth auxiliary variable v is always zero, 
the corresponding points are non-hyperbolic, and thus one should apply the center manifold 
analysis [122] in order to deduce whether it is unstable, saddle or stable. However, since 
all the above points have z c = 0, v is completely decoupled from the autonomous system 
(3.24)-(3.27). Therefore, one can examine only the three eigenvalues 1/1,1^2,1^3 displayed in 
Table 8, and the results of the stability analysis are presented in Table 4. Finally, in order 
to examine correctly the stability of points Fq and O2 we use the same linear approximation 
with that in the end of appendix A.l (for points Cq and 0\), and we deduce that Fq is a 
stable node while 0\ is saddle, with eigenvalues depicted in Table 8. 

We end this appendix with the following comment. The above stability conditions for 
F, Fq and G we extracted using the fact that the ^-direction is decoupled. However, as we 
mentioned the complete dynamical picture should be obtained after a full center manifold 
analysis [122], which is the correct way to handle zero eigenvalues. Therefore, in the fol- 
lowing we present this analysis for completeness. Since compact variables are preferable to 
see the subtleties on the dynamics, we perform the calculations using a Poincare projection 
method introducing the variables 

x r = pcos 6 sin <1> simp, y r = pcos ip, z r = psin 6 sin <I> sin^, v r = pcos $ sin^, (B.l) 



where p = ^7=5, r = \J x 1 + y 2 + z 2 + v 2 , [0,2vr], and -f < V < § [116-119]. 

• Center manifold analysis for F 

Using the Poincare coordinates (B.l), the point F transforms to the point with co- 

A 1 / 

ordinates x r = —-^j^iUr = 2 V ^ f ' Zr = ®i v r = 0. Introducing the new variables 



u = v r , 
vi = z r , 



v, = -^V 6 - ^ + » ( 2 - -f ) - for _ 2) - V 2 - f ■ 



1 r— r. I. \l\ Wi-#A 4 „ /r~A? 



2 



V 3 = \x r ^-\ 2 y\v + \yr (Ay - 6) + jZ r \j 1 " ~^~^V + \ \j 2 " ^ , (B.2) 

we deduce that the local invariant center manifold of F is given by the approximated 
graph 



(u,v 1 ,v 2 ,v 3 )\v 1 = (u 4 ) , V2 = - l -^2-^u 2 + (n 4 ) 



V3 = y2-^-u 2 + 0{u 4 ),\u\ <s\, (B.3) 
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where 5 > is a suitably small positive real number and 0(4) denotes terms of fourth 
order in the vector norm. Therefore, the dynamics at the center manifold is given by 



u = V2u 2 X v , \u\ < 5, (B.4) 
with solution for u(0) = no given by 

1 ^ 1 — y/2uQ Ayr ^ 

Since the origin is an inflection point of the potential U(u) = — ^-u 3 Xy of the 
gradient-like equation (B.4), it is implied that F is a saddle point for Xy < 3, only 
for perturbations in the direction of tv-axis and thus for perturbations in the u-axis. 



Center manifold analysis for F( 



o 



Using the Poincare coordinates (B.l), the point Fq transforms to the point with 
coordinates x r = 0, y r = z r = 0, v r = 0. Proceeding similarly to point F above 
and introducing the new variables 

ui = v r , vi = x r , v 2 = y r - -jp ^3 = z, (B.6) 
we find that the center manifold is given by the graph 

Uuuv u v 2 ,v 3 ) : Vl = -^= + 0(4),v 2 = 0(4), v 3 = 0(4), \ Ul \ < s\ , (B.7) 

where 5 is a suitably small constant and 0(4) denotes terms of fourth order in the 
vector norm. The dynamics on the center manifold is governed by the equation = 
0(4), from which it follows that Fq is indeed stable, but for n ^ not asymptotically. 
This behavior is depicted in Fig. 5. 

Center manifold analysis for G 

Using the Poincare coordinates (B.l), the point G transforms to the point with co- 
ordinates x r = t=£ y = z r = 0, v r = 0. Introducing the new variables 

u = v r , 

V\ = Zr, 



x r (2\ v -3\ 2 v + 9) \v{y/24-7\l + \v)+6 

V2 = ~ , — , 

y r (-At- + 12X 2 V + 3^24 - 7XyXv + y/24 - 7X 2 V X V + 9 
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Figure 5. We consider the coordinate system (B.6), using Ay = and n arbitrary (for the 
numerics we choose n = 1 but different n's correspond to the same projection) . The dashed thick 
line represents the center manifold of Fq, which is stable. The orbits above and below that line tend 
to it as time goes forward. 
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with (5 > a suitably small positive real number. The dynamics at the center 
manifold is given by 
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with solution for u(0) = uq given by 

«(r) = U ° X J (B.ll) 

X 2 v -3u ^X 2 v + 3t 

u 3 A2 _|_3 

Since the origin is an inflection point of the potential U(u) = — °f the 

gradient-like equation (B.10), it follows that G is a saddle point for either 3 < Xy < -y 
or Xy > only for perturbations in the direction of iv-axis and thus for perturba- 
tions in the v-axis. 

B.2 Stability of the critical points at infinity 

We consider the Poincare variables (B.l) (we restrict the angle ip to this interval since the 
physical region is given by y > 0), and thus the points at infinity (r — > +00) are those 
having p — > 1. Furthermore, the physical phase space is given by 

(x r ,y r ,z r ,v r ) £ [-1,1] x [0,1] x [-1,1] x [-1,1] 

such that 

x 2 + y 2 x 2 z r [\^6nv r x r — 6v 2 — 6x 2 — 6y 2 — 6z 2 + 6) 



l-v 2 -x 2 - y 2 - z 2 , (1 _ v 2 _ x 2 _ y 2_ z 2f/2 



< 1, (B.12) 



and v 2 + x 2 + y 2 + z 2 . < 1. 

Performing the transformation (B.l) on the system (3.24)-(3.27), and taking the limit 
p — > 1, the leading terms in the resulting system are 

(1 - p 2 )p —t --n(n + 1) cos 2 9 sin 2 $ cos 2 $ sin 4 tjj [2cos(2#) sin 2 $ sin 2 ip + cos 2^ + l] , 

(B.13) 

(1 - p 2 ) 2 9' -> -^n(n + l)sin6»cos 3 sin 2 (2$) sin 4 V>, (B.14) 

(1 - p 2 ) 2 & — > n{n + 1) cos 2 9 sin 2 $ cos 2 $ sin 5 ip cos ip [003(26*) sin 2 $ - l] , (B.15) 
(1 - p 2 ) 2 ^ 1 -> n(n + 1) cos 2 9 cos(2(9) sin $cos $sin tp. (B.16) 

Since the equation for p decouples, we only need to investigate the subsystem of the angular 
variables. Thus, equating (B.14), (B.15) and (B.16) to zero we obtain four classes of 
solutions: 
Class A: 

x r = 0, y r = costp, z r = ±sin$sinV ; , v r = cos $ sin ijj, (B.17) 

for some ip. These points correspond to 9 = ±S. For these points, £Ide of (3.36) diverges 
as p — > 1~, unless y r = 0, thus the meaningful critical points at infinity are those contained 
in the invariant circle z 2 + v 2 = 1. Representative critical points of this curve are the points 

L : [x r = 0, y r = 0, z r = ±1, v r = 0] 
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and 

M ± : [x r = 0, y r = 0, z r = 0, v r = ±1] . 

Class B: 

x r = 0, y r = costp, z r = 0, v r = ±sinip, (B.18) 

for some tp. These points correspond to $ = 0,7r, but since Qde of (3.36) diverges as 
p — > 1~ unless y r = 0, all these points are unphysical apart from 

M ± : [x r = 0, y r = 0, z r = 0, v r = ±1} . 

Class C: 

x r = 0, y r = 1, z r = 0, v r = 0. (B.19) 

These points correspond to ip = 0. They lead to a divergence in Qde m (3.36) and thus 
they are unphysical. 
Class D: 

x r = ± cos 9 sin?/;, y r = cos ip, z r = zbsin^sin^, v r = 0. (B.20) 

These points correspond to $ = but since &de of (3.36) diverges unless x r = y r = 0, 
all these points are unphysical apart from 

: [x r = 0, y r = 0, z r = ±1, v r = 0] . 

From the above analysis we deduce that for the singular points at infinity satisfy 

x r = (we mention that x r = is not an invariant set for the flow unless y r = or A = 0) 
that is points M^, or v r = that is points L^. 

Analyzing the dynamical system for x r ,y r , z r ,v r around M^- , we find that the eigen- 
values of the linearization are (0, |, — |, |), and thus, although non-hyperbolic, points 
behave as saddle points for the Poincare vector field. 

In order to calculate the eigenvalues of the linearization of the dynamical system for 
x r ,y r , z r ,v r around L^, we proceed as follows: having calculated the linear perturbation 
matrix Q we evaluate it at y r = v r = 0, then we take the limit x r — > and finally we 
incorporate the lateral limit z r — > V (that is from below) for L + or z r — > — 1 + (that is 
from above) for L~ . Calculating the eigenvalues of the resulting matrix we obtain that 
the eigenvalues of the linearization for both are given by (3, |, |). Thus, we deduce 
that points are unstable. 
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